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Abstract 

We show that the space of density matrices for ri-qubit 
states, considered as a (2"')^-diniensional real vector 
space, has a basis consisting of density matrices of sta- 
bihzer states. We describe an apphcation of this resuh 
to automated verification of quantum protocols. 

1 Definitions and Results 

We are working with the stabilizer formalism [S], in 
which certain quantum states on sets of qubits are repre- 
sented by the intersection of their stabilizer groups with 
the group generated by the Pauli operators. The sta- 
bilizer formalism is defined, explained and illustrated in 
a substantial literature; good introductions are given by 
Aaronson and Gottesman [1] and Nielsen and Chuang 
[a Sec. 10.5]. 

In this paper we only need to use the following facts 
about stabilizer states. 

1. The standard basis states are stabilizer states. 

2. The set of stabilizer states is closed under applica- 
tion of Hadamard (H), Pauli (X, Y, Z), controlled 

'l 0^ 



not (CNot), and phase (P 







) gates. 



Lemma 2 //O ^ a;, y < 2" and x ^ y then \x) + \y) and 
\x) + i\y) are stabilizer states. 

Proof: By induction on n. For the base case {n — 1), 
the closure properties imply that |0) -|- |1), |0) and 
1) + i\0) = |0) — i|l) are stabilizer states. 

For the inductive case, consider the binary represen- 
tations of X and y. If there is a bit position in which x 
and y have the same value b, then jcc) -|- |y) is the tensor 
product of \b) with an (n — l)-qubit state of the form 
\x') + |y'), where x' ^ y' . By the induction hypothesis, 
W) + \y') is a stabilizer state, and the conclusion follows 
from the closure properties. Similarly for \x) + i\y). 

Otherwise, the binary representations of x and y are 
complementary bit patterns. In this case, \x) + \y) can 
be obtained from GHZ„ by applying X to certain qubits. 
The conclusion follows from Lemma [T] and the closure 
properties. The same argument applies to \x) + i\y), 
using iGHZ„. □ 

Theorem 1 The space of density matrices for n-qubit 
states, considered as a {2'^)'^ -dimensional real vector 
space, has a basis consisting of density matrices of n- 
qubit stabilizer states. 

Proof: This is the space of Hermitian matrices and its 
obvious basis is the union of 



3. The set of stabilizer states is closed under tensor 
product. 

Notation 1 Write the standard basis for n-qubit states 
as {\x) I ^ X < 2"}, considering numbers to stand 
for their binary representations. We omit normalization 
factors when writing quantum states. 

Definition 1 Let GHZ„ = |0} + |2" - 1) and iGHZ„ 

|0) + i\2" — 1), as n-qubit states. 

Lemma 1 For all n, GHZ„ and iGHZ„ are stabilizer 
states. 

Proof: By induction on n. For the base case {n = 1), we 
have that |0) + |1) and |0) are stabilizer states, by 

applying H and then P to |0). 

For the inductive case, GHZ„ and iGHZ„ are obtained 
from GHZ„_i ® |0) and iGHZ„_i |0), respectively, by 
applying CNot to the two rightmost qubits. □ 



{|a;)(a;| | sC x < 2"} 



(1) 



{\x){y\ + \y){x\\0^x<y<2^} (2) 
{-i\x){y\+i\y){x\\0^x<y<2''}. (3) 
Now consider the union of 



{\x){x\ |0<x<2"} 



(4) 



{{\x) + \y))i{x\ + {y\)\0^x<y<2"} (5) 
{i\x) + i\yMx\ - i{y\) \0^x<y< 2"}. (6) 

This is also a set of (2")^ states, and it spans the space 
because we can obtain states of forms ([2]) and ([3]) by 
subtracting states of form from those of forms (O 
and (O. Therefore it is a basis, and by Lemma it 
consists of stabilizer states. □ 
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The stabihzer formahsm can be used to implement an 
efficient classical simulation of quantum computation, if 
quantum operations are restricted to those under which 
the set of stabilizer states is closed — i.e. the Clifford 
group operations. This result is the Gottesman-Knill 
Theorem [6]. Gay, Nagarajan and Papanikolaou [3j|4l|8] 
have applied it to formal verification of quantum sys- 
tems, adapting model-checking techniques 2^ from clas- 
sical computer science. A simple example of model- 
checking is the following. 

Consider a quantum teleportation protocol as a sys- 
tem with one qubit as input and one (different) qubit 
as output; call this system Teleport. Also consider the 
one-qubit identity operator /. Then the specification 
that teleportation should satisfy is that Teleport = I, 
where equality means the same transformation of a one- 
qubit state. The aim of model-checking in this context 
is to automatically verify that this specification is satis- 
fied, by simulating the operation of the two systems on 
all possible inputs. For this to be possible, the telepor- 
tation protocol is first expressed in a formal modelling 
language analogous to a programming language. 

The approach of Gay, Nagarajan and Papanikolaou re- 
duces the problem to that of simulating teleportation on 
stabilizer states as inputs, which can be done efficiently 
because the teleportation protocol itself only uses Clif- 
ford group operations. Correctness of teleportation on 
stabilizer states is interpreted as evidence for, although 
not proof of, correctness on arbitrary states. Now, how- 
ever, we can draw a stronger conclusion. 

The teleportation protocol defines a superoperator; 
this can be proved, for example, by the techniques of 
Selinger [S], who uses superoperators to define the se- 
mantics of a programming language that can certainly 
express teleportation. Superoperators, among other 
properties, are linear operators on the space of density 
matrices. To check equivalence of two superoperators, it 
is therefore sufficient to check that they have the same 
effect on all elements of a particular basis. By taking a 
basis that consists of stabilizer states, this can be done 
efficiently. 

Moreover, the number of stabilizer states on n qubits 
is approximately 2'" whereas the dimensionality of 
the space of n-qubit density matrices is only (2")^ = 2^". 
It is therefore more efficient to model-check on a basis 
than on all stabilizer states. 

3 Conclusion 

We have proved that the space of n-qubit density ma- 
trices has a basis consisting of stabilizer states, and ex- 
plained how this result can be used to improve the effi- 
ciency and strengthen the results of model-checking for 
quantum systems. 
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